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We report results of an acoustic duration reproduction task with stimulus duration of 2, 4, 
and 6 s, using 45 emotionally negative, positive, and neutral sounds from the International 
Affective Digitized Sounds System, in a sample of 31 young healthy participants. To 
investigate the influence of induced emotions on perceived duration, the effects of 
emotional modulation were quantified in two ways: (1) via model-free indices (aggregated 
ratios of reproduced times), and (2) via dual klepsydra model (DKM)-based estimates 
of parameters of internal time representation. Both data-analytic approaches reveal an 
effect of emotional valence/arousal, namely, a significantly longer reproduction response 
for emotional stimuli than for the neutral stimuli. The advantage of the DKM-based 
approach is its ability to disentangle stimulus-related effects, which are represented 
by "flow intensities," from general effects which are due to the lossy character of 
temporal integration. We explain the rationale of the DKM-based strategy and interpret 
the observed effect within the DKM-framework as transient increase of internal "flows." 
This interpretation is in line with recent conceptualizations of an "embodiment" of time 
where the model-posited flows correspond to the ongoing stream of interoceptive (bodily) 
neural signals. Neurophysiological findings on correlations between the processing of 
body signals and the perception of time provide cumulative evidence for this working 
hypothesis. 



Keywords: arousal, dual klepsydra model, duration reproduction, emotion, time perception 



1. INTRODUCTION 

Perception of duration is known to be dependent on many fac- 
tors, from physical characteristics of perceived events (Goldstone 
and Goldfarb, 1963; Block, 1978; Grondin, 1993) to psychophys- 
iological states of the perceiving subject (Wittmann, 2009; 
Mella, Conty, and Pouthas, 2011; Droit-Volet et al., 2013a). 
Experimental procedures and data-analytic models that allow us 
to disentangle and selectively test the role of different modulating 
factors are required for a better understanding of processes under- 
lying time perception, and temporal consciousness in general. 

In the present study we investigated the perception of tem- 
poral intervals marked by acoustic stimuli of varied emotional 
character, using the duration reproduction paradigm. The theo- 
retical framework of the study was the hypothesis of bodily states 
being the physiological basis for time perception (Wittmann, 
2009; Wittmann et al., 2010a). This hypothesis is based on the 
embodiment approach by Craig (2009) who proposes that our 
perception of time relates to emotional and visceral processes that 
all share a common underlying neural system within the inte- 
roceptive system and the insular cortex. In this context, several 
studies have shown how emotions and bodily arousal lead to an 
overestimation of duration. The experience of duration thereafter 
emerges from the processing and representation of emotional and 



body states (Droit-Volet and Gil, 2009; Droit-Volet et al., 2013b; 
Wittmann, 2013). This hypothesis fits conceptually with mod- 
els of internal representation of temporal duration based on the 
integration of intra-organismic "flows" which may be tentatively 
identified with streams of ascending neural signals (Meissner 
and Wittmann, 2011; Sysoeva et al, 2011). Specifically, we use a 
"klepsydraic" model, belonging to the class of "lossy integration" 
models, to quantify, test and interpret observed effects. 

The "dual klepsydra model" (dkm) (Wackermann and Ehm, 
2006) assumes that duration of an attended time interval is rep- 
resented by the state of an inflow-outflow unit (IOU), acting as a 
"lossy integrator" (see Appendix A for details). In the duration 
discrimination or reproduction task, two integrators are allo- 
cated, one being active ("filled") during the first interval ("encod- 
ing" phase), the other during the second interval (comparison or 
reproduction phase). Importantly, neither the "inflows" entering 
the integrators nor the states of the integrators are observable; 
they thus do not provide a direct measure of "subjective time." 
Only the result of their comparison is accessible to the sub- 
ject's awareness. In a duration reproduction task, equality of the 
states of the two integrators is indicated by the subject's verbal 
or motor response (e.g., key press) according to the instruction 
(Figure 1). The integrators are subject to continuous outflow, at 
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FIGURE 1 | Dual klepsydra model of duration reproduction, s = 

encoded duration, w = inter-stimulus interval, r = reproduced duration. 
[From Wackermann and Ehm (2006) with permission from Elsevier 
Science.] 



a rate proportional to their momentary states. Due to this "lossy" 
character of the lOUs, reproduced duration r is not a veridi- 
cal reproduction but rather a non-linear function of "encoded" 
duration s [Appendix A, Equation A3]. The DKM thus naturally 
accounts for the progressive shortening of the relative reproduc- 
tion response (i.e., decrease of the ratio r/s with increasing s), 
which is an ubiquitous phenomenon in duration reproduction 
experiments in the supra-second range (Spati, 2005; Wackermann 
and Ehm, 2006; Wackermann et al., 2008; Piitz et al., 2012). 

The DKM shares the concept of integrator (sometimes called 
"accumulator") with the "pacemaker-counter" model (PCM), 
which is influential in the time perception literature (Zakay and 
Block, 1997). In the DKM the input to integrators are continu- 
ous "flows," while the PCM postulates series of discrete "pulses," 
but this does not make a substantial difference: counting is only a 
special case of integration. The essential difference between the 
two models consists in their operating principles. In the DKM 
approach, different areas of the neural substrate may be allocated 
and re-used as integrators (cf. lech et al, 2005), their states are 
(bio)physical states, and their magnitudes compared directly. By 
contrast, the PCM is based on an idea of a central "counter," states 
of which are transferred to auxiliary "registers," and from there 
compared numerically with memory data. Briefly, the DKM is a 
sort of primitive, analog computer implemented in the biologi- 
cal matter, while the PCM mimics digital computing machinery 
(Wackermann, 2012). 

Yet another difference between the two models is that the DKM 
makes it possible to separate, conceptually as well as computa- 
tionally, effects due to the lossy character of neural representation 
(outflow) from effects caused by variation of the integrator input 
(inflow). The PCM, in its classic form, could account for effects 
of all kind only in terms of variation of the effective pulse 
frequency — i.e., number of "pulses" incoming at the accumulator 
per (physical) time unit (Treisman, 1963) — unless the memory or 
decision components are conceived of as active structures, such 
as in variants of the Scalar Timing models where the representa- 
tion of duration is modified (Gibbon, Church and Meek, 1984; 



Wearden, 2004). In the DKM, the same unit serves as integrator 
and memory storage, and the internal representation results from 
two concurrent processes, inflow and outflow. This is impor- 
tant in the context of the present study, where we aimed at 
changes in subjectively perceived duration induced by experimen- 
tal variations of the emotional valence of acoustic stimuli marking 
temporal intervals. To assess these effects, we have to disentangle 
them from the omnipresent "progressive shortening" effect men- 
tioned above. In the DKM this is possible, because the two kinds 
of effects are represented by two model parameters, k, regulating 
the outflow rate, and t), representing inflow ratios between the 
encoding and reproduction phase (Appendix A). 

The aim of the present study was to explore the effects of emo- 
tional modulation on duration reproduction in the supra-second 
range (2-6 s), using acoustic stimuli from a standardized, valence 
and arousal level rated system (IADS; Bradley and Lang, 2007). So 
far only a few studies have employed affective sounds in the anal- 
yses of temporal processes using the IADS (Noulhiane et al, 2007; 
Mella, Conty, and Pouthas, 2011). With this approach we aimed at 
extending evidence for the relationship between time perception 
and the processing of acoustic emotional content. In addition, we 
applied a data-analytic procedure based on the dual klepsydra 
model to differentiate valence/arousal-related effects from gen- 
eral, emotion-unspecific effects. With this strategy we present a 
novel parametric approach to analyze effects of varied stimulus 
quality on perceived duration. 

2. MATERIALS AND METHODS 

2.1. PARTICIPANTS 

Thirty-one subjects (16 women), students from the local univer- 
sity between 19 and 34 years of age (mean = 25.2 years; S.D. = 
3.2 years) were recruited for this study. All participants reported 
good health and no known neurological or psychiatric problems. 
They signed an informed consent before the beginning of the 
session and were reimbursed for taking part in the study. 

2.2. DESIGN AND STIMULI 

The International Affective Digitized Sounds System (IADS-2, 
Bradley and Lang, 2007) consists of 167 naturally occurring 
sounds that are rated according to the emotion dimensions of 
"valence," "arousal," and "dominance." Forty-five stimuli were 
selected from the IADS-2 and assigned to three groups of 15 
sounds each, according to their valence rating 

• Positive (p), e.g., pleasant music, a boy laughing, erotic sounds; 

• Negative (n), e.g., babies crying, a person vomiting, an alarm 
clock ringing; and 

• Neutral (e), e.g., clicking of type-writer, rain, animal noises. 

The positive (pleasant) and negative (unpleasant) sounds were 
selected such that they were comparable regarding average arousal 
rating. The reference numbers of the stimuli and statistics of 
ratings for the stimulus groups are shown in Table 1. 

This relatively large amount of stimuli was chosen in order 
not to present a given stimulus twice to a subject. The experi- 
ment consisted of 45 trials = the number of selected sounds ( 1 5 
of each group). Three different durations, 2, 4, and 6 s, were used 
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Table 1 | Acoustic stimuli used in the study and their average ratings. 



Group IADS no. Pleasantness Arousal 

Mean SD Mean SD 

POSITIVE 

a 200, 254, 802, 817 820 6.58 6.27 

b 110,351,366,717,810 7.22 5.78 

c 220,311,201,355,816 7.14 6.60 

Total 6.98 0.51 6.22 0.96 



NEGATIVE 

a 252,293,600,611,709 3.33 5.96 

b 255,288,296,712,730 2.82 6.35 

c 250,260,282,711,732 3.29 6.51 

Total 3.15 0.59 6.27 0.98 



NEUTRAL 

a 113, 246, 705, 720, 724 4.97 4.41 

b 102, 130, 376, 700, 701 4.87 4.55 

c 322,361,627,722,728 4.80 4.67 

Total 4.88 0.17 4.54 0.38 



in the duration reproduction task. This range of durations was 
chosen because the acoustical stimuli necessitate certain duration 
for being meaningful and evoking an emotion. Being too short 
they might not be deciphered correctly and when being too long, 
they might evoke modulations in subjective reactions. The stimuli 
have been validated for 5 s duration. 

Accordingly, each of the three stimulus groups (positive, neg- 
ative, neutral valence) was split into three sub-groups of five 
sounds each (a, b, c), which were assigned to the three duration 
conditions. The complete factorial structure of the experiment 
was thus 3 emotional valences (P, N, E) x 3 durations (s = 2, 4, 6) 
x 5 stimuli per sub-group. 

2.3. PROCEDURE 

The experiment was controlled by a portable computer, using 
a program written in Matlab (Mathworks Inc., Natick, MA) 
in conjunction with the Psychophysics Toolbox extensions for 
PC (Brainard, 1997). Stimuli were presented via headphones 
(Sennheiser HD 201). Stimuli were presented to all subjects with 
the same loudness, preset at a comfortable hearing level. 

Each trial started with the presentation of one of the IADS-2 
stimuli for 2, 4, or 6 s (encoding phase). After a pause of either 1.5, 
2, or 3 s, a pure sinewave tone of 150 Hz frequency was presented. 
Subjects had to stop the tone by pressing the space-bar when they 
felt that its duration had reached the duration of the first stim- 
ulus (reproduction phase). The next trial started after a variable 
inter-trial interval with durations of either 3, 3.5, or 4 s. Subjects 
received four practice trials before the proper experimental phase. 

Since it is tempting for subjects to use chronometric count- 
ing as a strategy to be more accurate, participants were strongly 
requested not to count and not to use any other strategies (e.g., 
subvocal singing) to measure the duration of sounds. It was 
emphasized that we were not interested in their counting abilities 



but in their subjective impression of duration. We did not employ 
a secondary task as it may interfere with the actual timing task. It 
has been shown that merely instructing subjects not to count is 
the best strategy for assessing interval timing abilities (Rattat and 
Droit- Volet, 2012). 

3. DATA REDUCTION AND ANALYSIS 

3.1. DATA PREPROCESSING 

Following the procedure described above, data from 45 trials were 
obtained from each subject, with five trials for each combination 
of three emotional valences (v = N, P, e) and three different pre- 
sented durations (s; = 2, 4, 6 s for i = 1, 2, 3, respectively). Time 
elapsed from the onset of the second tone to the space-bar press 
was measured as the subject's response r v ,u- These data were 
reduced to arithmetic means, 

! 5 

r v,t = - ^ r v,i,, (1) 

;=i 

for further analyses. (To simplify the notation, the superscripted 
bar will be omitted in the following text. The symbol r v j is to be 
understood as the average across a five-trial sub-group for given 
valence v and duration s,.) 

In the subsequent sections we apply two parallel strategies in 
order to aggregate data over different durations s,. The first strat- 
egy is purely descriptive, representing effects of varied stimulus 
valence by means of simple arithmetic indices; the second strat- 
egy is based on parametrization of the effects via the DKM. In this 
way we separate a straightforward demonstration of the existence 
of expected effects (model-free analysis) from their more refined 
quantification and interpretation (model-based). 

3.2. AGGREGATE RATIO INDICES 

To assess the effect of emotional valence v = N, P with respect to 
(w.r.t.) the neutral condition E, we calculated for each subject the 
mean ratio of reproduction responses 

taking the average over all three durations. If stimulus valence has 
no effect on the response, the expected value of the ratios r Vi ,-/r Eii -, 
and therefore of the aggregated ratios a v , are 1 (Appendix B). This 
will be tested across the sample of 31 subjects using a one-sample 
f-test 1 . 

3.3. DKM-BASED ANALYSIS 

To assess the effect of emotional valence v = N, P w.r.t. to the neu- 
tral condition E, we estimate inflow ratios between the encoding 



Obviously, this is a fair test as long as differences r„ ; — r E , are relatively 
small. However, for larger differences we may have a "built in" asymmetry in 
the evaluation procedure, since negative differences are mapped to the finite 
domain (0, 1) while positive differences are mapped to the infinite domain 
(i, oo). To symmetrize the distribution, we apply log-transformation to the 
a v values and run a t-test on log a v against the expected value of 0. 
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and reproduction phase, and from these we calculate net inflow 
ratios between the emotionally laden and neutral stimuli. (See 
Equations 6, 7 in Appendix C, where also a detailed rationale for 
the analysis is given). For each individual subject, the procedure 
consists of the following steps: 

1. Estimate k from the merged data-set consisting of all 45 trials, 
assuming formally r\ = 1. 

2. Using the estimate of k from step 1, estimate r\ v separately for 
v = P, N, E, each subset consisting of 15 trials. 

3. Using the estimates from step 2, calculate ratios r| PE = r|p /r) E 
andr) NE = WIe- 

Similar to the model-free analysis described in the preceding sub- 
section, the null-hypothesis-based expected value of inflow ratios 
T) PE and r] NE are 1. We proceed in the same vein, using one-sample 
f-tests to test for deviations from the expected value. 

4. RESULTS 

4.1. DESCRIPTIVE STATISTICS 

The grand means and standard deviations of reproduced dura- 
tion calculated across all 31 subjects are summarized in Table 2. 
The statistics of reproduced duration are shown as functions of 
the encoded duration s = 2, 4, and 6 s, sorted by emotional con- 
ditions P (positive), N (negative) and E (neutral), and for all trials 
regardless of the emotional condition (the bottom row). 

As a global test of significance a Two-Way ANOVA was com- 
puted, with factor 1 = Valence (three levels) and factor 2 = 
Duration (three levels). The results (Greenhouse-Geisser cor- 
rected) are: 

• Factor Valence: F = 8.089 (df = 2, 60); P < 0.001. 

• Factor Duration: F = 308.075 (df = 2, 60), P < 0.0001 (triv- 
ial). 

• Interaction Valence versus Duration: F = 3.809 (df = 4, 120); 
P = 0.012. 

Calculating post hoc f-tests for dependent variables (df = 30) 
within duration conditions and between emotional conditions, 
we found 



• for s = 6 s, significant differences between conditions N and E 
(t = 2.986, P < 0.01), as well as between conditions P and E 
(f = 2.226, P < 0.05). 

In addition, a review of the reproduction responses relatively to 
the encoded duration s reveals, on the average, a minor (non- 
significant) over-reproduction for s = 2 s, and a pronounced 
under-reproduction for s = 4 and 6 s (Figure 2). We thus observe 
the progressive shortening effect, consistent with other duration 
reproduction studies (cf. Introduction). 

4.2. AGGREGATE RATIO INDICES 

Empirical distributions of a N , a P in the sample of 31 subjects are 
shown in Figure 3; sample means are 1.0175 and 1.0715, respec- 
tively. One-sample f-test against the expected value of 1 yields 
f = 1.171 (ns) for the condition N, and t = 4.646 (P < 0.001) for 
the condition P (df = 30 for both tests). We thus have a highly 
significant lengthening effect of ~7% for positive valence stim- 
uli; in addition, we observe a small (1.75%) but not significant 
difference in the same direction for negative valence stimuli 2 . 

4.3. DKM-BASED ANALYSIS 

Results of k estimates based on KRF fits to merged data-sets were a 
necessary pre-requisite for further analysis; therefore these results 
are reported first. In four cases we obtained an estimate of k = 0. 
Data inspection revealed that these four participants consistently 
and sometimes vigorously over-reproduced the presented dura- 
tion s. This kind of response, sometimes seen in a minority of 
subjects (cf. Ptitz et al, 2012), cannot be accounted for by the 
DKM. A negative value of k would be required to fit these data, but 
k < 0 is biophysically meaningless; in such cases the estimation 
procedure stops at the lower-bound k = 0. 

Furthermore, for one subject we obtained an estimate of k = 
0.24 s _1 . Visual inspection of her data revealed that her responses 
r varied in a very narrow range and only minimally reflected the 
s durations presented in the encoding phase. It is unclear whether 
she misunderstood the instruction or exhibited an unusual kind 



2 The f-tests applied to log-transformed indices yield t = 0.951 (ns) for the 
condition N, and t = 4.582 (P < 0.001) for the condition P, thus confirming 
the findings based on a v indices. 



• For s = 4 s, a significant difference between conditions P and E 
(f = 2.600, P < 0.05); 



Table 2 | Statistics of reproduced durations (grand means and 
standard deviations) sorted by emotional valence and encoded 
duration s. 



Valence 


s = 


2s 


s = 


4s 


s 


6s 


Mean 


SD 


Mean 


SD 


Me; m 


SD 


Positive 


2.335 


0.690 


3.717 


0.869 


4.799 


1.102 


Negative 


2.136 


0.571 


3.456 


0.876 


4.839 


1.106 


Neutral 


2.210 


0.676 


3.496 


0.937 


4.527 


1.212 


All 


2.227 


0.646 


3.556 


0.892 


4.722 


1.137 



A 

r[s] 



i mean, all conditions 




r[s] 
6 



6 s[s] 



o positive 
• negative 
o neutral 




6 S[S] 



FIGURE 2 | Average reproduction response for the entire sample of 
31 subjects. (A) All conditions (shown are grand means ±1SD). (B) 
Separated by conditions E,N,P (SD bars omitted for the sake of legibility). 
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of cognitive bias. In any case, this extreme value, largely deviating 
from the rest of the sample, is to be considered an evident outlier. 

The four over-reproducing subjects and the one non- 
responder were removed from the data-base; the cleaned sample 
thus consists of 26 subjects. In this sample, k values were in the 
range from 0.003 to 0.096 s _1 , median k = 0.03 s _1 . 

Empirical distributions of the inflow ratios T] NE and T] PE in 
the sample of 26 subjects are shown in Figure 4; sample means 
are 1.0421 and 1.0788, respectively. One-sample f-test against 
the expected value of 1 yields t = lAii (P < 0.05) for r| NE , and 
t = 4.081 (P < 0.001) for n PE (df= 25 for both tests). We thus 
observe increase of inflow ratios w.r.t. the neutral stimuli in both 
emotion conditions: a moderate (~4%) and significant effect for 
negative valence stimuli, and a large (~8%) and highly significant 
effect for positive valence stimuli 3 . 



3 The f-tests applied to log-transformed inflow ratios yield ( = 1.905 (P ~ 
0.06) for fine ' and * = 3 - 958 ( p < 0.001) for r| PE . In spite of somewhat 
reduced significances, the results are in agreement with those obtained for 
original inflow ratios. 



5. DISCUSSION AND CONCLUSION 

Our results complement earlier studies showing that emotional 
stimuli are judged to last longer than more neutral stimuli of the 
same physical duration (Droit-Volet and Gil, 2009; Wittmann, 
2009). An overestimation of duration is more likely seen with 
negative emotional valence. However, the presentation of posi- 
tive emotional stimuli can also lead to a relative overestimation 
of duration as compared to neutral stimuli (Lambrechts et al., 
201 1), results which our data supports. In fact, in our task with 
the specific selection of sounds from the IADS battery, positive 
sounds are relatively longer reproduced than negative sounds, 
both emotional stimuli being relatively over-estimated compared 
to neutral sounds. Such findings of an overestimation are typically 
discussed within the framework of an increased physiological 
arousal level that leads to a higher frequency of an internal pace- 
maker underlying time perception (Wittmann and Paulus, 2008; 
Gil and Droit-Volet, 2012). A higher speed of an internal clock 
would lead to the accumulation of more pulses emitted by a 
hypothetical pacemaker operating in the seconds range (Zakay 
and Block, 1997). For example, in one study a faster tempo in 
music, which was related to an increased subjective arousal, pro- 
duced relative duration overestimates as compared to slow tempi 
(Droit-Volet et al., 2013b). In another study, temporal intervals 
were judged to last longer after subjects had viewed a frightening 
film that increased fear and arousal levels (Droit-Volet, Fayolle 
and Gil, 201 1). In this line of evidence, looming stimuli that were 
virtually approaching the viewer were relatively overestimated as 
compared to steady stimuli or those that were virtually receding 
(Wittmann et al, 2010b). 

In the present analytical approach we want to frame the results 
related to the temporal overestimation of emotionally arousing 
stimuli within a novel framework provided by the dual klepsy- 
dra model of duration representation. Wackermann and Ehm 
(2006, p. 489) suggested that external stimulus-induced, phasic 
effects are mediated by changes in integrator inflows (assessed by 
parameter r\), whereas effects of the global organismic state and 
its slow, tonic changes, are mediated by the outflow rate (assessed 
by parameter k). There is now substantial empirical support for 
this conjecture. Parameter k has been shown to be a test-retest sta- 
ble, probably genetically co-determined, individual characteristic 
(Sysoeva et al., 2010). This does not mean that k is perfectly con- 
stant; it is subject to intra- individual circadian variations (Spati, 
2005), and can be affected by neurochemical agents acting on the 
synaptic level (Wackermann et al., 2008). All these findings point 
out to the slow, tonic character of k reactivity, on time scales of 
hours or longer. On the other hand, the inflow ratio r\ has been 
shown to reflect physical properties of stimuli used in a time per- 
ception task (brightness-duration interaction; Wackermann and 
Blankenburg, 2009). It is thus reasonable to expect that the inflow 
ratio r) is the proper "locus of effect" also for stimulus-induced 
influences of other kinds. 

The working hypothesis of the present study is based on the 
conjecture of two different loci of the effects discussed above. 
Stimuli of high emotional valence — and, usually of high biologi- 
cal relevance — induce a general physiological response (arousal), 
resulting in a temporary increase of interoceptive afferences. We 
hypothesize that the observed effects on time perception are due 
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FIGURE 3 | Distributions of aggregated ratio indices a P and a N in the 
entire sample (n = 31 subjects). Upper part: individual data points; below: 
estimates of probability density functions (thin curve: a^, thick curve: ap), 
obtained with a Gaussian kernel (Silverman, 1986) of bandwith o = 0.04. 
Vertical dotted line indicates the null-hypothesis-based expectation. 



o-o o-oo-keccoo — cxyo O-OD— 



— o-o-o— co-o-coo— o©-cco<D-ooo- 




1-5 n V E 



FIGURE 4 | Distributions of estimated flow ratios n PE and ii NE in the 
cleaned sample (n = 26 subjects). Upper part: individual data points; 
below: estimates of probability density functions (thin curve: T)p E , thick 
curve: tine), obtained with the same method as in Figure 3. Vertical dotted 
line indicates the null-hypothesis-based expectation. 
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to a temporary increase of neural flows feeding the hypotheti- 
cal integrators. Therefore, phasic changes of arousal level should 
be measurable in terms of inflow ratios between arousing versus 
emotionally neutral conditions. To obtain robust and valid esti- 
mates of these changes, we aggregate reproduction data over the 
entire range of encoded duration s. This is an important differ- 
ence from the majority of duration reproduction studies testing 
effects for different values of s separately. In our approach it is the 
response function r = f (s) in its entirety that carries information 
on the internal duration representation. Accordingly, instead of 
seeking/testing effects in the data space, we first apply a mathe- 
matical model to the data and then test for the expected effects in 
the parameter space. 

In the present paper we used two data-evaluation strate- 
gies in parallel, to prevent a possible criticism for moulding 
the data to our preferred model. The results of both strategies 
are well comparable: Pearson's correlation between a N and r| NE 
is R = +0.857 (df= 24, P < 0.001), and between a ? and r| PE 
it is R = +0.878 (df= 24, P < 0.001). We can thus say that 
the aggregated ratio indices (model-free) and the inflow ratios 
(model-based) measure almost the same — yet not exactly the 
same thing. The DKM-based strategy turns out to be more sensi- 
tive, particularly w.r.t. effects induced by negative -valence stimuli, 
seemingly non-significant when using the ratio indices, but sig- 
nificant when quantified in terms of inflow ratios 4 . Arguably, this 
better sensitivity is due to separation of stimulus-specific effects 
from the unspecific effect of "progressive shortening." The latter 
effect is demonstrated by the summary statistics of reproduction 
responses (Table 2), showing an average shortening of ~1 1 % for 
s = 4s, and ~21 % for s = 6s. This effect is further evidenced 
by the distribution of individual k values which are positive for 
the great majority of subjects; the sample median value of k = 
0.03 s _1 is in very good agreement with k values found in earlier 
studies (Spati, 2005; Wackermann and Ehm, 2006; Wackermann 
and Spati, 2006; Wackermann et al., 2008; Sysoeva et al., 2010; 
Wittmann et al., 2011; Piitz et al, 2012). 

The main finding of the present study is a lengthening of 
reproduction response for emotional stimuli w.r.t. to the response 
for neutral stimuli: we found significant effects for both condi- 
tions P and N versus E. The two dimensions are not independent 
in our selection of stimuli: both stimulus groups P and N show 
on the average higher arousal than the group E (cf. Table 1). 
Therefore we cannot clearly differentiate between influences of 
the two stimulus dimensions; nonetheless, our results are in line 
with the general arousal-increase hypothesis specified above. 

The main methodological innovation of the present study is 
the use of inflow ratios to quantify the stimulus-induced effect on 
duration perception. We wish to emphasize that our data-analytic 



4 To substantiate this claim, we calculated statistics of the aggregated ratios 
after removing the five subjects we had excluded from the DKM analysis: mean 
ii N = 1.0169, mean flp = 1.0758. One-sample t-tests against 1 yield t = 1.024 
(ns) for the condition N and t = 4.720 (P < 0.001) for the condition N (df = 
25 for both tests), differing only minimally from those obtained for the whole 
sample of 31 subjects (section 4.2). Since these results are based on the same 
subset of 26 subjects that entered the DKM analysis, we conclude that the better 
performance of the DKM-based strategy is not an artifact due to the difference 
between the two samples. 



approach makes these effects not only detectable — i.e., testable 
for their presence and statistical significance — but also measur- 
able and comparable on a ratio scale. The quantification can work 
both ways: the method allows for characterization of stimuli of 
different kinds by their potence to alter perceived duration in 
a representative sample of observers. Or, inversely, susceptibility 
of individual observers to emotion-induced effects can be mea- 
sured, using a standardized stimulus set, as in our study. Our 
data do not allow us to assess the intra-individual stability of 
those measures; repeated sessions and test-retest analyses would 
be required for that purpose. However, the notion of susceptibil- 
ity to emotion-induced alteration of time experience, considered 
as individual trait, is supported by high intra-individual correla- 
tions between measures obtained separately for P and N stimuli: 
Pearson's correlation between a N and a P is -R = +0.566 (df = 
29, P < 0.01, see Figure 5), and notably higher, R = +0.744, 
between r| NE and r| PE (df= 24, P < 0.001, see Figure 6). This 
higher intra-individual correlation obtained with the DKM-based 
parameters is an additional argument in favor of our model-based 
approach 5 . 

Another interesting, even if only qualitative, observation con- 
cerns the distribution of r| NE ratios (Figure 4). In fact, the local 
maximum at n NE ~ 1.1 almost exactly coincides with that for 
r|p E , and the sample mean of rj NE = 1.04 thus somewhat under- 
estimates the effect of negative-valence stimuli. This is due to 
a secondary protrusion in the r| NE distribution, seen at values 
slightly below 1. The shape of the distribution thus suggests a 
mixture of two unimodal distributions, that of "responders" of 



5 After removing the five subjects that were excluded from the DKM analysis, 
the correlation between «n and a P becomes slightly insignificantly higher: R = 
+0.589 (d/= 24,P < 0.01), still clearly below the intra-individual correlation 
between the DKM-based flow-ratio estimates. 
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FIGURE 5 | Correlation plot of aggregated ratio indices a N versus a P 
(n = 31, /?= +0.566). 
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(n= 26, R= +0.744). 







the same average level as to the positive-valence stimuli, and that 
of "non-responders" to negative-valence stimuli. A deconvolution 
of the r| NE distribution would be interesting, but such analysis is 
not possible because of a rather small sample size. 

Concluding, the present study (1) confirmed and extended 
evidence for emotion-induced alterations of duration percep- 
tion, (2) demonstrated the practical utility of the dual klepsydra 
model for quantification of such effects, and (3) corroborated 
the hypothesis of internal time representation as resulting from 
integration of flows of interoceptive signals reflecting emotionally 
modulated bodily states. Saying that our approach opens a "win- 
dow on" neurophysiology of subjective time might be unjustly 
exaggerated; integration of behavioral response measures with 
simultaneous physiological measurements would be desirable. It 
is not exaggerated to say, however, that the quantitative meth- 
ods designed for the purpose of the present study open at least a 
"peephole" into the integration mechanisms underlying duration 
representation, and subjective time awareness in general. 

ACKNOWLEDGMENTS 

Marc Wittmann was supported by the European project COST 
ISCH Action TD0904 "Time In MEntaL activitY: theoreti- 
cal, behavioral, bioimaging and clinical perspectives" (TIMELY; 
www.timely-cost.eu). We wish to thank Oksana Gutina for con- 
ducting part of the experiments and for general assistance. We are 
thankful to two reviewers of the original manuscript for critical 
remarks and constructive suggestions. 

REFERENCES 

Block, R. A. (1978). Remembered duration: effects of event and sequence complex- 
ity. Mem. Cogn. 6, 320-326. doi: 10.3758/BF03197462 

Bradley, M. M., and Lang, R J. (2007). The International Affective Digitized Sounds 
(2nd edn; IADS-2): Affective Ratings of Sounds and Instruction Manual (Technical 
report B-3). Gainesville, FL: University of Florida. 



Brainard, D.H. (1997). The psychophysics toolbox. Spat Vis. 10, 433-436. doi: 

10.1163/156856897X00357 
Craig, A. D. (2009). Emotional moments across time: a possible neural basis for 

time perception in the anterior insula. Phil. Trans. R. Soc. B 364, 1933-1942. 

doi: 10.1098/rstb.2009.0008 
Droit- Volet, S., Fayolle, S. L., and Gil, S. (2011). Emotional state and time 

perception: effects of film-induced mood. Front. Integr. Neurosci. 5:33. doi: 

10.3389/fnint.2011.00033 
Droit- Volet, S., Fayolle, S., Lamotte, M., and Gil, S. (2013a). Time, emotion and the 

embodiment of timing. Timing Time Percept. 1, 99-126. doi: 10.1163/22134468- 

00002004 

Droit-Volet, S., and Gil, S. (2009). The time-emotion paradox. Philos. Trans. R. Soc. 
B 364, 1943-1954. doi: 10.1098/rstb.2009.0013 

Droit-Volet, S., Ramos, D., Bueno, J. L. O., and Bigand, E. (2013b). Music, emotion 
and time perception: the influence of subjective emotional valence and arousal? 
Front. Psychol. 4:417. doi: 10.3389/fpsyg.2013.00417 

Gibbon, J., Church, R. M., and Meek, W. H. (1984). "Scalar timing in memory," in 
Timing and Time Perception, eds J. Gibbon and L. Allan (New York, NY: New 
York Academy of Sciences), Ann. N.Y. Acad. Sci. 423, 52-77. 

Gil, S., and Droit-Volet, S. (2012). Emotional time distortions: the fundamental 
role of arousal. Cogn. Fmot. 26, 847-862. doi: 10.1080/02699931.2011.625401 

Goldstone, S., and Goldfarb, J. L. (1963). Judgment of filled and unfilled 
durations: intersensory factors. Percept. Mot. Skills 17, 763-774. doi: 
10.2466/pms.l963. 17.3.763 

Grondin, S. (1993). Duration discrimination of empty and filled intervals 
marked by auditory and visual signals. Percept. Psychophys. 54, 383-394. doi: 
10.3758/BF03205274 

Jech, R., Dusek, P., Wackermann, J., and Vymazal, J. (2005). Cumulative 
blood oxygenation level- dependent signal changes support the 
'time accumulator' hypothesis. NeuroReport 16, 1467-1471. doi: 
10.1097/01.wnr.0000175616.00936.1c 

Lambrechts, A., Mella, N., Pouthas, V, and Noulhiane, M. (2011). Subjectivity 
of time perception: a visual emotional orchestration. Front. Integral Neurosci. 
5:73. doi: 10.3389/fnint.2011.00073 

Meissner, K. and Wittmann, M. (2011). Body signals, cardiac aware- 
ness, and the perception of time. Biol. Psychol. 86, 289-297. doi: 
10.1016/j.biopsycho.201 1.01.001 

Mella, N., Conty, L., and Pouthas, V. (2011). The role of physiological arousal 
in time perception: psychophysiological evidence from an emotion regulation 
paradigm. Brain Cogn. 75, 182-187. doi: 10.1016/j.bandc.2010.11.012 

Noulhiane, M., Mella, N., Samson, S., Ragot, R., and Pouthas, V. (2007). How 
emotional auditory stimuli modulate time perception. Emotion 7, 697-704. doi: 
10.1037/1528-3542.7.4.697 

Ptitz, P., Wittmann, M., and Wackermann, J. (2012). Duration reproduc- 
tion: lossy integration and effects of sensory modalities, cognitive func- 
tioning, age, and sex. Percept. Motor Skills Learn. Mem. 115, 370-384. doi: 
1 0.2466/22. 10.27.PMS.1 15.5.370-384 

Rattat, A. C, and Droit-Volet, S. (2012). What is the best and easiest method of 
preventing counting in different temporal tasks? Behav. Res. Meth. 44, 67-80. 
doi: 10.3758/sl3428-01 1-0135-3 

Silverman, B. W. (1986). Density Estimation for Statistics and Data Analysis. 
Monographs on statistics and applied probability, Vol. 26. London: Chapman 
& Hall. 

Spati, J. P. (2005). Time Perception in the Context of Circadian Activity. Diploma 
thesis, Department of Neurobiology and Biophysics, University of Freiburg, doi: 
10.1007/978-1-4899-3324-9 

Sysoeva, O. V., Tonevitsky, A., and Wackermann, J. (2010). Genetic determinants of 
time perception mediated by the serotonergic system. PLoS ONE 5:el2650. doi: 
10.1371/journal.pone.0012650 

Sysoeva, O. V., Wittmann, M., and Wackermann, J. (2011). Neural representation 
of temporal duration: coherent findings obtained with the 'lossy integration' 
model. Front. Integr. Neurosci. 5:37. doi: 10.3389/fnint.201 1.00037 

Treisman, M. (1963). Temporal discrimination and the difference interval: impli- 
cations for a model of the 'internal clock'. Psychol. Monogr. 77, 1-31. doi: 
10.1037/h0093864 

Wackermann, J. (2012). "On clocks, models and metaphors: understanding the 
klepsydra model," in Multidisciplinary Aspects of Time and Time Perception, eds 
A. Vatakis, A. Esposito, M. Giagkou, F. Cummins, and G. Papadelis (Berlin: 
Springer), 246-257. 



Frontiers in Neurorobotics 



www.frontiersin.org 



February 2014 | Volume 8 | Article 11 | 7 



Wackermann et al. 



Effects of emotion on duration reproduction 



Wackermann, J., and Ehm, W. (2006). The dual klepsydra model of internal 
time representation and time reproduction. /. Theor. Biol. 239, 482-493. doi: 
10.1016/j.jtbi.2005.08.024 

Wackermann, J., and Meyer-Blankenburg, N. (2009). "Differential effect of stim- 
ulus brightness in visual duration discrimination," in Fechner Day 2009, eds 
M. A. Elliott, S. Antonijevic, S. Berthaud, P. Mulcahy, B. Bargary, C. Martyn, 
and H. Schmidt (Galway: International Society for Psychophysics), 459-464. 

Wackermann, J., and Spati, J. (2006). Asymmetry of the discrimination function 
for temporal durations in human subjects. Acta Neurobiol. Exp. (Wars.) 66, 
245-254. 

Wackermann, J., Wittmann, M., Hasler, E, and Vollenweider, F. X. (2008). Effects 
of varied doses of psilocybin on time interval reproduction in human subjects. 
Neurosti. Lett. 435, 51-55. doi: 10.1016/j.neulet.2008.02.006 

Wearden, J. (2004). Decision processes in models of timing. Acta Neurobiol Exp. 
(Wars.) 64, 303-317. 

Wittmann, M. (2009). The inner sense of time. Philos. Trans. R. Soc. Lond. B 364, 
1955-1967. doi: 10.1098/rstb.2009.0003 

Wittmann, M. (2013). The inner sense of time: how the brain creates a representa- 
tion of duration. Nature Rev. Neurosci. 14, 217-223. doi: 10.1038/nrn3452 

Wittmann, M., and Paulus, M. P. (2008). Decision making, impulsivity, and time 
perception. Trends Cogn. Sci. 12, 7-12. doi: 10.1016/j.tics.2007.10.004 

Wittmann, M., Simmons, A. N., Aron, J. L., and Paulus, M. P. 
(2010a). Accumulation of neural activity in the posterior insula 
encodes the passage of time. Neuropsychologia 48, 3110-3120. doi: 
10.1016/j.neuropsychologia.2010.06.023 



Wittmann, M., van Wassenhove, V., Craig, A. D., and Paulus, M. P. (2010b). The 

neural substrates of subjective time dilation. Front. Human Neurosci. 4:2. doi: 

10.3389/neuro.09.002.2010 
Wittmann, M., Simmons, A. N., Flagan, T., Lane, S. D., Wackermann, J., and Paulus, 

M. P. (2011). Neural substrates of time perception and impulsivity. Brain Res. 

1406, 43-58. doi: 10.1016/j.brainres.2011.06.048 
Zakay, D, and Block, R. A. (1997). Temporal cognition. Curr. Direct. Psychol Sci. 6, 

12-16. doi: 10.1111/1467-8721.epll512604 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 26 November 2013; accepted: 11 February 2014; published online: 28 
February 2014. 

Citation: Wackermann J, Meissner K, Tankersley D and Wittmann M (2014) Effects 
of emotional valence and arousal on acoustic duration reproduction assessed via the 
"dual klepsydra model". Front. Neurorobot. 8:11. doi: 10.3389, fnbot.2014.00011 
This article was submitted to the journal Frontiers in Neurorobotics. 
Copyright © 2014 Wackermann, Meissner, Tankersley and Wittmann. This is an 
open-access article distributed under the terms of the Creative Commons Attribution 
License (CCBY). The use, distribution or reproduction in other forums is permitted, 
provided the original author(s) or licensor are credited and that the original publica- 
tion in this journal is cited, in accordance with accepted academic practice. No use, 
distribution or reproduction is permitted which does not comply with these terms. 



Frontiers in Neurorobotics 



www.frontiersin.org 



February 2014 | Volume 8 | Article 11 | 8 



Wackermann et al. 



Effects of emotion on duration reproduction 



r[s] 




0 10 20 30 s[s] 

r[s] 




FIGURE A1 | Upper part: different forms of the KRF for n = 1 (equal 
inflows) and k varied from 0 to 0.04 with step 0.01. The higher the loss 
coefficient k, the more pronounced is the KRF curvature, i.e., the 
progressive shortening with increasing duration s. Lower part: different 
forms of the krf for a constant value of k = 0.03 and r| varied from 0.8 to 
1.2 with step 0.1. The curvature is basically determined by k and further 
modulated by the inflow ratio r). — ISI in both displays are constant 0.5 s. 
Thin diagonal line indicates the chronometrically correct response r = s. 



APPENDIX 

A DUAL KLEPSYDRA MODEL: THE ESSENTIALS 

A linear inflow-outflow unit (iou) is described by an ordinary 
differential equation 

dy 

-f=i- K y, (Al) 
df 

where y denotes the momentary state, i is the "inflow" inten- 
sity, and k > 0 is a constant. For i being a piecewise constant 
function of time t, the IOU state y is a piecewise exponential 
function. The dual klepsydra model (dkm), designed specifi- 
cally for the duration reproduction paradigm, consists of two 
such lOUs (Figure 1). During the encoding phase (0 < t < s), 
IOUi is filled with intensity i\ > 0, and i\ = 0 afterwards. During 
the reproduction phase (s + w < f), IOU2 is filled with intensity 
z'2 > 0, until the equality of states, y\ (f) = /2(f) is reached and the 
subject's response elicited at time t = s + w + r. The following 
relation holds for durations r, s and w. 



Denoting, for short, r\ = and rewriting Equation (A2) to 
express the reproduced duration r as a function of the encoded 
duration s and the pause w between the two intervals, we obtain 

r = K- 1 ln(l + Ti (l-e^ 5 ) e- KW ) , (A3) 

known as the "klepsydraic reproduction function" (krf). Its form 
is determined by two system parameters, k and t), and an exper- 
imental parameter w. The KRF is a monotonically increasing, 
negatively accelerated function of s, approaching asymptotically 
a limiting value for s — > 00. Parameter k measures the intrinsic 
loss of the hypothetical neural integrators; its physical dimension 
is [k] = s . The inverse value k -1 is proportional to the half- 
outflow time of the IOU with zero inflow. The larger is k (i.e., the 
shorter the characteristic time k -1 ), the stronger is the negative 
curvature of the KRF, and the more pronounced is the progressive 
shortening effect in the reproduction data (Figure Al). 

In experiments employing stimuli of the same physical or 
psychological characteristics for the encoding and reproduction 
phase, T) = 1 is usually assumed. A KRF fit to reproduction data, 
using e.g., the weighted least squares method (Wackermann and 
Ehm, 2006), yields an estimate of the parameter k. Reciprocally, 
if k is known, and experimental conditions are varied between 
the encoding and reproduction phase, the corresponding inflow 
ratio r| between conditions can be calculated (Wackermann and 
Blankenburg, 2009). Simultaneous fitting of k and r| to a sin- 
gle data-set is not recommendable: estimation errors for k and 
t), respectively, are strongly correlated (Wackermann and Ehm, 
2006) and simultaneous estimates of the two parameters are thus 
unreliable. 

For k — > 0 (loss-less integration) the KRF becomes a linear 
function of s, r = r\ s, that is, a deviation from veridical repro- 
duction can only be assessed in terms of flow ratios. In that case 
the DKM is functionally equivalent to the PCM. 



B RATIONALE FOR THE AGGREGATE RATIOS ANALYSIS 

Tests for effects of emotional valence of presented stimu- 
lus w.r.t. the neutral condition E can be based on dif- 
ferences or ratios between average responses. For our pur- 
pose the use of ratios r v /r E is preferable, as these ratios 
are naturally linked to the "pacemaker-counter" model (PCM) 
(Zakay and Block, 1997). 

Let/i iV denote the pulse frequency during encoding of the first 
interval, s, in a condition v. The ratio between the reproduced 
duration r and the encoded duration s then gives an estimate 
of the inverse ratio between the respective frequencies, r v /s = 
f\,v/fi- The second stimulus used in the reproduction phase being 
always the same, we assume a constant pulse frequency fa in both 
conditions. 
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Thus we have 



r v /s 
r E /s 



/l,E 



(Bl) 



This means that the between-conditions ratio of the reproduced 
durations provides a direct estimate of the between-conditions 
ratio of the pulse frequencies. The aggregate indices a v (Equation 
2) used in our analyses provide estimates of the respective fre- 
quency ratios averaged across all durations s,- (i = 1, 2, 3). 

C RATIONALE FOR THE DKM-BASED ANALYSIS 

Our analysis was based on the hypothesis that emotional valence 
of presented acoustic stimuli modulates the inflows entering the 
neural integrators. Emotional valence was manipulated during 
the encoding phase; assume that for an individual subject i\ 
attains three different values, i ? , z N , or i E , depending on the stim- 
ulus. The second stimulus used in the reproduction phase was 
always the same, a sine-wave reference tone, so we assume that the 
inflow intensity ii is constantly i K {. The inflows per se are inob- 
servable; only data-based inferences on inflow ratios are possible. 
For the three emotion conditions we define 



_ If _ ! n 

t]p = - — > tin = - — ■ 

Zref hd 



TIE = 



!ref 



(CI) 



Of interest now are effects of emotionally laden stimuli, P or N, 
w.r.t. the neutral stimulus E. For this purpose we take ratios 



_ Tip _ h _ t)n _ ! N 

T)pe = — — , T| NE = — — 

T)e l E T) E Z E 



(C2) 



providing net inflow ratios between emotion versus neutral con- 
ditions. 

A prior knowledge of individual k is needed to estimate T| v for 
any of the three conditions v = P, N, E. Ideally, k should be esti- 
mated from a data-set with perfectly homogeneous stimuli (e.g., 
reference tone used for both encoding and reproduction), but 
our experimental design does not yield such data. As an heuris- 
tic approximation we estimate k from a merged data-set, i.e., not 
distinguishing between emotion conditions and putting formally 
t) = 1. This is not quite correct, because 2/3 of the data were 
obtained for emotion-versus-reference stimuli, and so the subse- 
quent estimates of r\ v (Equation CI) may be somewhat biased. 
However, we are mainly interested in T) vE ratios: as seen from 
Equation (C2), the unknown intensity i K ( cancels out, and it is 
thus legitimate to expect that the procedural bias is compensated 
by this operation. 

For k — > 0 (loss-less integration) the ratio r v /s gives a direct 
estimate of r\ v , by analogy to the PCM (cf. Equations Bl, C2). 
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